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Foreword 


Numerous  naval  applications  are  directly  dependent  on  knowing  the  past, 
present,  and  future  sea  surface  temperature  (SST)  fields.  The  Diurnal  Ocean 
Surface  Layer  (DOSL)  model  at  the  Fleet  Numerical  Oceanography  Center 
attempts  to  fill  a  portion  of  these  requirements  by  forecasting  the  24-hour 
change  in  global  SSTs.  This  procedure  is  done  by  using  the  wind  stress  and 
heat  flux  forcing  functions  available  from  the  Naval  Operational  Global 
Atmospheric  Prediction  System  as  input  to  the  parameterized  DOSL  model. 

This  report  deals  with  validating  the  DOSL  products  using  infrared  imagery 
obtained  in  several  locations  during  the  night  and  day  to  map  the  diurnal  SST 
changes.  Particular  attention  ic  focused  on  extreme  diurnal  warming  events, 
since  these  cases  are  most  commonly  associated  with  surface  mixed-layer 
warming  that  disrupts  the  surface  acoustic  duct.  Several  events  are  detailed 
in  the  subtropical  Atlantic  to  depict  the  capabilities  within  the  DOSL  model. 


Tpi'Jly' 

W.  B.  Moseley 
Technical  Director 


J.  Tupaz,  Captain,  USN 
Commanding  Officer 


Executive  Summary 


The  Diurnal  Ocean  Surface  Layer  (DOSL)  model  at  the  Fleet  Numerical 
Oceanography  Center  forecasts  the  24-hour  change  in  global  sea  surface 
temperatures  (SST).  Validating  the  DOSL  model  is  a  difficult  task  due  to  the 
huge  areas  involved  and  the  lack  of  in  situ  measurements.  Therefore,  this  report 
details  the  use  of  satellite  infrared  multichannel  SST  imagery  to  provide  day 
and  night  SSTs  that  can  be  directly  compared  to  DOSL  products.  This  water- 
vapor-corrected  imagery  has  the  adva'itages  of  high  thermal  sensitivity 
(0. 12°C),  large  synoptic  coverage  (nearly  3(X)0  km  across),  and  high  spatial 
resolution  that  enables  diurnal  heating  events  to  be  readily  located  and  mapped. 

Several  case  studies  in  the  subtropical  North  Atlantic  readily  show  that  DOSL 
results  during  extreme  heating  periods  agree  very  well  with  satellite-imagery- 
derived  values  in  terms  of  the  pattern  of  diurnal  warming.  The  low  wind  and 
cloud-free  conditions  necessary  for  these  events  to  occur  lend  themselves  well 
to  observation  via  infrared  imagery.  Thus,  the  normally  cloud-limited  aspects 
of  satellite  imagery  do  not  come  into  play  for  these  particular  environmental 
conditions. 

The  fact  that  the  DOSL  model  does  well  in  extreme  events  is  beneficial  from 
the  standpoint  that  these  cases  can  be  associated  with  the  destruction  of  the 
surface  acoustic  duct.  This  so-called  “afternoon  effect”  happens  as  the  after¬ 
noon  warming  of  the  mixed  layer  disrupts  the  sound  channel  and  the 
propagation  of  acoustic  energy.  Thus,  the  DOSL  model  can  be  used  to  help 
estimate  the  location  of  where  the  afternoon  effect  may  occur,  as  well  as 
help  analysts  interpret  artificially  induced,  diurnal  heating  SST  gradients  when 
generating  mesoscale  front  and  eddy  maps. 
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Diurnal  Ocean  Surface  Layer  Model  Validation 


I.  Introduction 

The  U.S.  Navy  requires  a  wide  range  of  real-time 
sea  surface  temperature  (SST)  analyses  (nowcasts)  and 
forecasts  to  meet  the  increasingly  sophisticated  needs 
demanded  by  sensors  and  weapon  systems.  The  SSTs 
may  directly  or  indirectly  impact  successful  and 
efficient  completion  of  duties  ranging  from  antisub¬ 
marine  warfare  (ASW),  optimum  track  ship  routing, 
and  search  and  rescue,  or  serve  as  vital  input  to 
atmospheric  prediction  models  that  require  accurate 
air-sea  interface  initial  conditions.  This  report  focuses 
on  the  validation  of  one  model  aimed  at  satisfying  the 
need  to  forecast  daily  SST  changes. 

The  Diurnal  Ocean  Surface  Layer  (DOST)  model 
(Clancy,  1986)  is  run  daily  at  the  Fleet  Numerical 
Oceanography  Center  (FNOC)  in  Monterey,  CA,  and 
predicts  the  daily  SST  changes  due  to  the  impact  of 
surface  heat  fluxes  and  wind  stress.  This  parameterized 
model  requires  no  SST  initial  conditions  and  uses  the 
forecast  winds  and  fluxes  from  the  Navy’s  global 
atmospheric  model  to  provide  the  key  input 
parameters.  It  then  produces  a  prediction  for  the 
diurnal  SST  heating  that  is  valid  globally  at  the  time 
of  local  maximum  heating. 

The  DOST  model  thus  generates  a  contoured  map 
for  both  hemispheres  outlining  the  regions  of  large 
diurnal  heating.  This  graphical  guidance  can  be  used 
to  outline  areas  where  the  surface  acoustic  duct  may 
be  lost  due  to  solar-heating-induced  stratification 
(afternoon  effect);  to  help  correctly  interpret  infrared 
SST  or  forward-looking  infrared  images  for  ocean 
mesoscale  feature  locations;  and  to  delineate  regions 
where  day  and  night  satellite  SSTs  should  not  be  mixed, 
since  they  represent  two  totally  different  environmental 
conditions.  The  DOST  model  can  thus  provide  some 
unique  capabilities  when  properly  used  in  conjunction 
with  the  right  Navy  mission. 

Validating  a  prognostic  SST  model  over  large  time 
and  space  scales  (DOST  uses  a  polar  stereo  grid  of 
~  320  km  and  is  interpolated  to  a  2.5°  x  2.5°  grid) 
raises  several  interesting  problems.  It  is  not  sufficient 
to  match  individual  point  values  with  in  situ  measure¬ 
ments,  since  the  pattern  of  SST  warming  is  also  a  key 
output  value.  Therefore,  accurate  SST  reports  within 
regional  and  sometimes  basin  scale  domains  must  be 


combined  with  the  attribute  of  coincident  timing 
needed  at  the  minimum  and  maximum  periods  of  solar 
heating.  Multichannel  SSTs  (MCSST)  from  polar 
orbiter  infrared  sensors  are  the  only  means  pr^^scntiy 
available  to  attempt  this  arduous  task. 

II.  Diurnal  Ocean  Surface  Layer  Model 

The  DOST  model  uses  the  sc,  'ing  relationships 
derived  by  Price  et  al.  (1986)  to  .elate  the  expected 
afternoon  mixed-layer  response  to  the  local 
atmospheric  forcing  A  model  similar  to  that  of  Price 
et  al.  (1978)  was  used  to  simulate  the  extensive  ocean 
mixed-layer  data  set  obtained  from  the  R/V  Flip  using 
rapid-profiling  conductivity-temperature-depth  and 
VMCM  current  meters.  These  detailed  observations 
and  model  runs  enabled  Price  et  al.  (1986)  to 
develop  and  verify  scaling  relationships  that  predict 
the  afternoon  SST  warming  and  transient  layer  depth 
using  surface  heat  fluxes  and  surface  wind  stress  values. 
This  parameterization  was  later  verified  by  Bowers 
et  al.  (1986)  and  Price  et  al.  (1987)  with  independent 
data  from  the  northwestern  Atlantic. 

These  scaling  relationships  are  the  basis  for  the 
DOST  model  now  functioning  in  real  time  with  global 
coverage  at  FNOC  (Clancy,  1986).  It  forecasts  the 
change  in  day-night  SST  (AT^  and  diurnal  thermocline 
depth  (D),  valid  everywhere  in  the  local  afternoon.  This 
diurnal  SST  variation  will  occur  locally  from  a 
minimum  near  0600  to  a  maximum  about  1400  local 
time.  The  change  due  to  diurnal  heating  is  given  by 
Equation  1: 

AT  =  1.5  ,  (1) 

C  S 

where  g  is  the  acceleration  of  gravity,  a  is  the  expansion 
coefficient  for  seawater,  C  is  the  heat  capacity  of 
.seawater,  H  is  the  maximum  net  downward  surface 
heat  flux  (latent  -i-  sensible  +  infrared  -i-  solar)  during 
the  day,  P  is  1/2  the  time  interval  over  which  the  net 
surface  heat  flux  is  downward,  S  is  the  wind  stress 
parameter,  and  F  is  a  radiation  extinction  function. 


The  parameter  S  is  calculated  from 

r«  _ _  /y’r»\i0.5 


/ 


[1  -  COS  (fP)] 


(2) 


where  t  is  the  surface  wind  stress  averaged  over  the 
period  of  net  downward  surface  heat  flux  and  /  is 
the  Coriolis  parameter.  The  radiation  function  F 
is  derived  from 


F  =  1 


0.4  (2/  -  H)  exp  i-kD) 


(3) 


where  /  is  the  maximum  surface  solar  heat  flux  during 
the  day  and  k  is  the  extinction  coefficient  for  solar 
radiation  (Clancy,  1986). 

The  forcing  functions  are  derived  from  the  Navy 
Operational  Global  Atmospheric  Prediction  System 
(NOGAPS)  model  (Rosmond,  1981).  This  model 
provides  forcing  functions  on  a  global,  73  x  144 
spherical  grid  that  yields  a  resolution  of  2.5°.  These 
forcing  fields  are  subsequently  interpolated  to  the 
63  X  63  polar  stereographic  grid  on  which  DOST  runs. 
This  latter  grid  provides  resolutions  that  vary  from 
200  km  (tropics)  to  nearly  400  km  (polar  regions).  In 
determining  S,  the  wind-speed-dependent  drag 
coefficient  of  Garrai  (1977)  is  used  to  calculate  surface 
wind  stresses  from  the  NOGAPS  19.5-m  winds.  Also, 
the  algorithms  of  Clancy  and  Poliak  (1983)  are  used 
to  interpolate  the  surface  heat  fluxes  and  wind  stresses 
between  the  6-hourly  values  provided  by  NOGAPS 
(Clancy,  1986). 

Particular  attention  should  be  paid  to  the  forcing 
function  grid  resolution.  The  subtropical  zones  with 
the  best  chances  of  strong  diurnal  heating  occur 
between  20°  and  30°  from  the  equator.  NOGAPS  has 
a  resolution  near  275  km  and,  thus,  cannot  clearly 
define  mesoscale  heating  events  that  may  occur  and 
cannot  be  verified  with  high-resolution,  satellite-derived 
SSTs.  The  forcing  function  spatial  definition  is,  thus, 
one  limiting  factor  in  the  DOST  accuracy. 

The  DOST  model  produces  a  forecast  once  per  day 
at  OOZ  (Greenwich  Mean  Time),  which  is  valid 
worldwide  in  the  local  afternoon.  This  means  that  the 
valid  DTG  (Date-Time-Group)  (86081500  is  OOZ 
15  Aug  86)  is  a  function  of  longitude.  The  valid  time 
at  a  particular  location  is 


Time  =  DTG  +  14  hr  -i- 


L  (Ihr) 
15° 


for  0  <  L  <  270°W 


(4) 


and 


Time  =  DTG  +  8  hr  + 


(Z,  -  270) 


for  270  <  /.  360°  W  , 


(.‘5) 


where  L  is  the  longitude  of  the  location  ranging 
from  0°  to  360°  measured  positive  westward  (W)  from 
Greenwich. 

The  table  illustrates  sample  valid  forecast  times  at 
different  longitudes  for  the  same  DTG  of  86081500. 


Longitude 

L  (°W) 

Valid  DTG  (Z) 

Local  Time 

15°W 

15° 

86081515 

1400 

75°W 

75° 

86081519 

1400 

150“W 

150° 

86081600 

1400 

150°  E 

210° 

86081604 

1400 

60°  E 

300° 

86081510 

1400 

15°E 

345° 

86081513 

1400 

Sample  graphical  DOST  output  (6  August  1987)  is 
shown  in  Figure  1.  The  large  FNOC  chart  size  and 
day/night  SST  difference  contours  at  0.2°C  intervals 
necessitate  exhibiting  just  the  Atlantic  and  East  Pacific 
basins:  otherwise,  interpreting  values  would  be  quite 
difficult.  One  can  readily  locate  areas  where  significant 
warming  (greater  than  1°C)  is  forecast  by  the  closely 
bunched  contours  in  numerous  regions. 

The  North  Atlantic  is  particularly  active,  with  greater 
than  1 .5°C  warming  predicted  by  the  DOST  model  for 
the  Gulf  of  Mexico,  the  Gulf  Stream,  the  Northeast 
Atlantic  south  of  Iceland,  and  the  two  local  maxima 
at  30°N  off  Northwest  Africa.  The  number  and 
intensity  of  changes  forecast  for  6  August  1987  is  not 
the  norm.  A  scan  through  weeks  of  charts  reveals 
that  the  patterns  change  dramatically  in  size,  in  shape, 
and  in  value  from  one  day  to  the  next  as  the  synoptic 
weather  pattern  evolves. 

Figure  2  illustrates  the  corresponding  surface  wind 
and  atmospheric  pressure  fields  (6  August  1986) 
available  from  the  FNOC  marine  wind  analysis.  The 
1024-mb  subtropical  high  centered  between  20°-30°N 
stretches  across  a  huge  longitudinal  range  and  readily 
dominates  the  weather  and,  thus,  cloud  cover.  This 
event  is  a  critical  factor  in  that  subsidence  within  the 
high  will  tend  to  deter  cloud  development,  thus  opening 
the  way  for  high  insolation.  Light  winds  are  also  the 
rule,  as  evidenced  by  the  lack  of  wind  barbs  within 
the  center.  (Higher  resolution  wind  fields  will  be 
covered  later.) 


III.  Diurnal  Warming  and 
SST  Measurements 

Insolation  affects  the  entire  mixed-layer  water 
column  to  various  degrees,  depending  on  several 
factors:  amount  of  solar  heating,  surface  wind  stress, 
water  clarity,  and  temperature  stratification  when  the 
heating  event  first  begins.  The  warming  is  initially  more 
abrupt  at  the  sea  surface  and  is  thus  reflected  in  the 
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Figure  I.  Sample  graphical  output  from  the  DOSL  model  for  6  August  1987.  Contours  every  0.2°C  represent  forecast  of  diurnal  change  in  SST 
valid  everywhere  at  local  time  of  maximum  insolation.  (Note  all  chart  values  should  be  multiplied  by  O.I.) 
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SST  as  a  “hot  spot”  when  compared  lo  surrounding 
ambient  SST  values  not  impacted  by  the  above  condi¬ 
tions.  As  these  favorable  diurnal  warming  conditions 
continue,  the  mi.xed  layer  shallows  and  warms,  creating 
a  stratified  layer  that  may  destroy  the  surface  acoustic 
duet . 

This  so-called  afternoon  effect  is  critical  to  Navy 
.•\S\\  strategics  (LTiek.  1983,  and  Laevastu  and 
Hubert.  1970).  Knowledge  of  when  the  diurnal 
warnring  w  ill  negate  use  of  the  surface  acoustic  duct  is 
important  for  a  variety  of  sensor/weapon  applications. 
The  DOSL  model  could  thus  be  used  as  one  of  several 
tools  (along  with  the  Thermodynamic  Ocean  Prediction 
■System — TOPS)  to  help  locate  where  the  afternoon 
effect  may  take  place.  .Appropriate  action  by  naval 
forces  could  then  lessen  the  chances  that  incorrect 
sensors  or  sensor  deployment  patterns  were  used  within 
the  specific  domain  of  interest. 

Intense  diurnal  warming  may  create  sharp,  hori¬ 
zontal  SST  gradients  that  could  be  readily  confused 
with  other  features  related  to  real  mesoscale  ocean 
frontal  boundaries.  The  generation  of  these  gradients  is 
especially  disconcerting,  since  the  regional  FNOC  three- 
dimensional  thermal  structure  analyses  (Optimum 
Thermal  Interpolation  System — OTIS)  are  gradually 
becoming  more  dependent  on  frontal  “bogu.ses” 
derived  by  analysts  who  extract  frontal  locations  and 
strengths  from  infrared  imagery.  Thus,  misleading  SST 
gradient  information  could  be  passed  onto  the  OTIS 
analysis  if  diurnal  warming  events  are  not  screened  out. 

Observations  of  sea  surface  warming  due  to 
insolation  (diurnal  heating)  have  been  inadequate  for 
many  years  (Cornillon  and  Stramma,  1985).  Ship 
weather  renorts  include  SST  readings  that  are  often 
representative  of  the  engine  intake  temperatures  (Earle, 
1985).  As  such,  they  occur  well  below  the  water  line 
(sometimes  as  deep  as  15-20  m)  and,  thus,  do  not 
indicate  the  true  SST.  Some  researchers  have  argued 
that  liiis  aspect  is  a  positive  factor  and  surface  skin 
effects  are  negated,  and  that  the  ship  values  are  a 
measure  of  the  “mixed-layer”  temperature. 

This  logic  does  not  hold  true  when  dealing  with  the 
topic  of  diurnal  surface  heating.  Surface  heat  fluxes 
and  wind  stress  are  the  key  to  determining  whether 
large  temperature  changes  will  occur  at  the  ocean 
surface  on  a  daily  time  scale.  Therefore,  the  skin 
temperature  (the  apparent  temperature  of  the  sea  that 
is  in  contact  with  the  air)  is  the  critical  sea  measure¬ 
ment,  not  some  temperature  value  at  a  depth  of  5  m 
or  10  m.  This  fact  significantly  raises  the  value  of 
satellite  infrared  SSTs  for  the  study  discussed  here. 

IV.  Saiellite-Derived  Sea  Surface 
Temperatures 

Satellite-derived  SSTs  have  progressed  significantly 
since  the  first  infrared  sensors  were  used.  Meteoro¬ 
logical  applications  paved  the  way  for  the  first  series  of 


environmental  satellites,  but  oceanography  was  a 
distant  by-product  due  in  part  to  poor  thermal 
sensitivity.  This  restriction  was  lifted  when  the 
TIROS-N  scries  of  National  Oceanic  and  .Atmospheric 
Administration  (NOAA)  polar  oribiters  began  using 
10-bit  digitization  (0.12  C  sensitivity)  over  the  whole 
range  of  possible  SSTs.  Small  SST  gradients  could  then 
be  viewed  and  surface  currents  mapped  as  never  before 
pos;:iblc.  However,  the  use  of  single-channel  infrared 
data  is  limited  by  the  attenuating  effect  of  atmospheric 
water  vapor.  Maul  (1983)  clearly  states  the  debilitating 
effect  water  vapor  can  have  by  calculating  potential 
SST  corrections  ranging  from  nearly  0  to  as  much  as 
lO^C.  These  values  are  truly  disconcerting  when 
realizing  that  such  errors  can  immediately  rule  out  the 
vast  majority  of  studies  that  rely  on  absolute  SST 
values  tor  quantitative  studies. 

In  1981,  NOAA  began  the  operational  production 
of  MCSSTs  by  using  the  differential  absorption 
observed  in  two  infrared  channels.  The  10.3-1 1.3)um 
and  11.5-12.5)im  channels  were  used  to  generate  a 
correction  value  applied  to  the  1 1 .5-12.5)im  value  that 
removed  the  major  water  vapor  effects.  Thus,  the 
differing  response  between  the  two  infrared  spectral 
bands  to  the  .same  columnar  water  vapor  permitted  the 
calculation  of  a  correction  factor,  which  fluctuates  as 
the  water  vapor  varies. 

The  NOAA  Advanced  Very  High  Resolution 
Radiometer  (AVHRR)  is  now  a  five-channel  sensor 
that  includes  an  additional  spectral  band  between 
3.5),4m  and  3.93)im.  This  addition  enables  the  use  of 
a  three-channel  algorithm  that  can  be  superior  to  the 
two-channel,  split-window  equations  (channels  4  and  5) 
when  the  ever-present  noise  in  channel  3  is  low  just 
after  launch  or  satellite  outgassing.  This  problem  will 
be  corrected  for  the  NOAA  K,L,M  series  by 
incorporating  slightly  different  sensor  technology. 

The  MCSST  allows  the  adjustment  to  be  done  on 
a  pixel-by-pixel  basis  and  eliminates  the  need  for 
additional  data  sources,  such  as  atmospheric  sounders, 
which  arc  plagued  with  poor  spatial  resolution  or 
radiosondes  that  provide  sparse  coverage  over  ocean 
domains.  NOAA  processes  global  area  coverage 
(GAC)  4-km  data  in  2  x  2  unit  arrays  to  produce 
8-km  MCSSTs  on  a  daily  basis.  The  procedure  to  create 
MCSSTs  has  been  refined  over  the  years,  and  their 
accuracy  has  continually  increased  to  the  point  where 
they  now  have  an  RMS  error  of  0.6-0.7°C  (Strong  and 
McClain,  1985)  while  producing  7(),(X)0  to  100,000 
global  retrievals  per  day.  This  measurement  milestone 
allows  MCSST  usage  for  a  wide  variety  of  research 
and  operational  projects. 

These  National  Environmental  Satellite,  i)ata. 
Information  Service  (NESDIS)  MCSSTs  arc  proc¬ 
essed  orbit  by  orbit  as  the  G  AC  data  are  downloaded 
from  the  spacecrtift  tape  recorders.  The  time  delay  from 
sensor  data  collection  to  MCSSf  retrieval  residing  in 


the  N[;S[)IS  cur'  day  ohscrvatii)n  tile  varies  from 
3  to  f)  hours.  I  hese  data  can  then  be  accessed  by 
operati(r.ial  activities  and  incorporated  into  SST  or 
three-dimensional  ocean  thermal  analyses. 

This  processing  procedure  is  designed  to  create  a 
reasonable  global  distribution  to  satisfy  a  variety  of 
SST  mapping  requirements.  However,  not  all  cloud- 
free  pixels  are  used  to  generate  MCSSTs.  The  MC'SST 
algorithm  incorporates  a  search  pattern  that  finds 

I  (low),  5  (medium)  or  15  (high)  observations  per 

II  »  II  CiAC  array  (about  50  x  50  km).  This  factor 
limits  the  total  number  of  MCSSTs  produced  and 
results  mainly  from  several  considerations — disk  space, 
processing  time,  and  overall  data  needs.  Thus, 
hundreds  of  thousands  of  potential  MCSST  retrievals 
are  not  produced  each  day. 

The  DOSl.  validation  effort  could  use  isolated 
matchups  between  day  night  MCSST  measurements 
and  DCSI.  model  output  but  is  unrealistic  for  several 
reasons:  no  DOST  digital  values  were  retained,  only 
the  contoured  charts;  interpolating  coarse  DOSL  values 
to  .MCSST  locations  may  induce  large  errors;  and  this 
method  would  not  permit  us  to  see  whether  the  pattern 
of  diurnal  heating  was  being  correctly  reproduced. 
Therefore,  .MCSSTs  (imagery)  provides  us  with  the 
synoptic  scale  coverage  and  highly  accurate  SST  values 
required  for  DOST  comparisons. 

V.  Validation  Technique 

Cloud  contamination  severely  restricts  any  validation 
effort  that  relies  solely  on  infrared  imagery.  Therefore, 
instead  of  acquiring  images  in  real  time,  the  project 
relied  on  the  DOSL  model  to  indicate  which  regions 
should  be  viewed  more  closely.  The  first  phase  revolved 
around  the  collection  of  TNOC’s  DOSl.  and  marine 
surface  wind  charts  for  the  1986  Northern  Hemisphere 
summer  months.  The  charts  were  subjectively  scanned 
to  select  only  events  that  exhibited  large  diurnal 
warming  over  a  spatial  domain  greater  than  lO'"  longi¬ 
tude  (approximately  KKK)  km). 

DOSL  charts  for  August,  September,  and  most  of 
October  1986  were  scrutinized  to  find  likely  maximum 
diurnal  heating  events  that  would  probably  be 
associated  with  few  clouds  and  thus  indicate  very  good 
infrared  imagery  v  iewing  conditions.  Dozens  of  cases 
were  located  in  the  North  Atlantic,  the  Mediterranean 
Sea,  and  the  central  and  eastern  North  Pacific,  as  well 
as  a  few  isolated  occurrences  in  other  spots.  These 
potential  case  studies  were  then  screened  by  taking  into 
account  the  conditions  usually  found  in  each  area. 

The  North  Pacific  is  well  known  for  the  extensive 
cloud  coverage  generated  by  the  combined  effects  of 
the  meteorological  and  oceanographic  conditions  that 
arc  manifested  through  air-sea  interaction  processes. 
These  conditions  are  prevalent  in  the  summer  and  the 
late  fall  when  visible  or  infrared  viewing  is  extremely 


poor  at  best.  However,  the  DO.SI  model  picJicled 
large  diurnal  warming  based  on  the  toieeast  ol  light 
winds  and  little  surface  mixing  even  though  skies  were 
cloudy.  The  North  Pacific  was  eliminated  from  further 
study,  since  the  infrared  data  were  not  simable. 

The  Mediterranean  Sea  was  initially  thought  to  be 
a  superior  location  for  testing  the  DOST  model,  but 
several  problems  arose.  Excellent  satellite  viewing  is 
the  rule  rather  than  the  exception  and  permits  regular 
infrared  surveys.  Light  winds  and  insolation  can  quite 
frequently  ctiuse  ihe  model  to  forecast  diurnal  heating 
of  2 'C.  However,  the  regions  of  heating  can  be  small, 
arc  quick  to  dissipate  and  reappear,  and  are  difficult 
to  time  with  regards  to  satellite  overpasses.  The  region 
also  frequently  contains  dust  and  aerosols,  which  can 
attenuate  infrared  signals  and  decrease  the  accuracy 
(vf  measurements. 

Earlier  efforts  by  Cornillon  and  Stramma  ( 1985)  and 
Stramma  et  al.  (1986)  showed  that  selected  regions  in 
the  Sargasso  Sea  experienced  diurnal  warming  in  excess 
of  1  "C  up  to  30  percent  of  the  time.  DOSL  charts  were 
checked  to  locate  days  when  significant  heating  was 
forecast  and  then  nearly  100  NOAA  AVHRR  digital 
images  vvere  scrutinized  to  find  cloud-free,  day/night 
image  pairs.  These  data  were  previously  collected  in 
real  time  by  the  NORDA  Satellite  Digital  Receiving 
and  Processing  System  (Hawkins  et  al.,  1985)  to 
support  I  he  GEOS  AT  Ocean  Application  Program 
(GOAP;  Lybanon  et  al.,  1987).  GOAP  combined  the 
synoptic  infrared  coverage  with  precision  GEOSAT 
altimeter  sea  surface  height  data  to  map  mesoscale 
front  and  eddy  features  within  the  Gulf  Stream  system. 
All  data  were  high  picture  resolution  transmission 
data  (1  km)  from  cither  Wallops  Island,  VA,  or  the 
NORDA  tracking  antenna. 

The  DOSL  model  output  indicated  numerous  cases 
where  warming  larger  than  DC  should  have  been 
visible  in  satellite  imagery.  However,  a  quick  look  at 
the  digital  images  revealed  very  few  suitable  cloud-free 
pairs.  .Many  day/night  MCSSl  difference  images  were 
generated,  but  none  were  of  sufficient  quality  toi  this 
study,  even  though  cloud-free  bits  and  pieces  did  verify 
well  with  DOSL  values. 

Attention  was  then  turned  to  the  central  and  eastern 
North  Atlantic,  since  the  Azores-Bermuda  high  was 
often  creating  good  cloud-free  conditions.  These 
conditions  were  evident  in  the  persistently  high  DOSL 
diurnal  heating  forecasts.  Numerous  AVHRR  digital 
images  were  ordered  from  NESDIS  for  further 
investigation.  Both  1-km  local  area  coverage  and  4-km 
GAC'  data  were  included  to  see  whether  or  not  one  was 
preferable  to  the  other  (note  that  Cornillon  and 
Stramma,  1985,  found  4  km  sufficient  for  their  study). 
Imagery  for  16-18  and  28  August;  3,  15,  23,  and 
30  September;  and  2  and  7  October  1986  were  then 
read  in  and  initially  scanned  for  cl'uid  contamination. 
A  quick  survey  concluded  the  August  data  were  far 
superior,  and  they  were  then  processed  further. 


The  small  number  ol  cases  were  deemed  insut'l'ieient 
m  quantity,  so  attention  was  turned  to  other  methods 
to  increase  the  ease  study  comparisons.  Lessons  learned 
during  the  1986  summer  were  then  incorporated  into 
a  near-real-time  et'tort  in  the  1987  summer.  DOSl. 
eharls  were  cheeked  daily  in  the  North  Atlantic  to  try 
and  catch  the  beginning  of  diurnal  heating  events. 
C  onditions  normally  changed  in  logical  phases  as  cloud 
systems  with  higher  winds  moved  out  and  were 
followed  bv  clearer  skies  and  weaker  winds.  Direct 
communication  between  FNOC  and  NORDA 
personnel  permitted  the  collection  of  several  more 
eases,  liaeh  example,  vsheiher  for  1986  or  1987,  was 
processed  in  the  following  manner; 

•  NO.'\A-9  channels  4  and  5  data  were  calibrated 
and  earth-located  using  software  on  NORDA’s  Intcr- 
aeiise  Digital  Satellite  Image  Processing  System 
(IDSIPS;  Hawkins  et  al.,  1985). 

•  liarth  location  was  fine-tuned  using  European, 
African,  and  island  landmarks  and  resulted  in 
registration  accuracies  of  2  pixels  or  better. 

•  k  hannels  4  and  5  were  used  for  both  day  and  night 
passes  instead  of  .L4,5  for  the  night  and  4,5  for  the 
day  to  be  consistent  and  to  reduce  the  introduction  of 
noise  inherent  in  channel  3.  Conservative  cloud  tests 
were  incorporated  that  likely  erred  on  the  side  of 
eliminating  some  cloud-free  pixels  to  make  sure 
partially  cloudy  pixels  were  not  allowed  to  slip  through 
and  introduce  erroneous  data  points  in  the  final  results. 

•  Multichannel  split-window  SSTs  were  calculated 
for  all  cloud-free  pixels  in  both  day  and  night  imagery 
using  the  NESDIS  operational  equations  valid  for  the 
period  (Strong  and  McClain,  1985). 

•  Daytime  MCSST  imagery  was  then  subtracted 
from  the  nighttime  versions  to  produce  a  day/night 
difference  image.  These  imagery  were  then  falsely 
colored  in  0.5'’C  increments  to  view  more  easily  the 
regions  ot  maximum  diurnal  heating.  Clouds  were 
masked  in  black  or  dark  gray. 

The  DOSl.  charts  were  not,  however,  in  an  easy 
format  to  compare  directly  with  the  images.  Thus,  the 
charts  were  digitized  on  IDSIPS  and  processed  to 
the  exact  image  map  projection  and  scale.  Isotherm 
contours  at  0.2°,  0.5°,  1.0°,  1 .5°,  2.0°,  and  2.5°C  were 
kept  and  all  others  were  interactively  erased  to  reduce 
the  cluttered  pattern  in  intense  heating  regions  (DOSL 
output  originally  contains  contours  at  0.2°C  intervals). 
DOSL  contours  could  then  be  directly  overlaid  in  a 
graphics  plane  on  top  of  the  falsely  colored,  day/nighl 
difference  images. 

VI.  Validation  Results 

August  1986  was  an  excellent  diurnal  heating  event 
time  due  to  the  meteorological  conditions  that  existed 
in  the  central  and  eastern  North  Atlantic.  The  sub¬ 
tropical  high  dominated  the  central  North  Atlantic  and 


contained  vasi  ar.'as  of  subsidence,  as  eudeiiced  m 
pressure  ma[is.  such  as  shovMi  in  1  igine  2.  Inirared 
imagery  was  relatively  cloud-free  for  a  large  region 
bounded  by  2.5-40°N  and  10-45  W  . 

A.  16  .Aujju.sl  1986 

f  igure  3  (1537Z  16  .August  1986)  n-  donnna'.cd  b> 
a  warm,  dark-red  patch  of  SS  I  s  oriented  nmiheast  to 
southwest  off  the  extreme  northwestern  Afiican  coast. 
This  section  of  very  warm  SS  fs  extends  westward  over 
20  longitude  and  is  cloudless  in  the  niosi  intense  aiea 
near  31°N  and  16  This  verv  warm  water  is  in  direst 
contrast  to  the  much  coolei  coastal  iipwclling  along 
the  northeastern  African  coast  and  contributes  to  a 
substantial  SST  gradient. 

The  day 'night  MCSST  image  tfig.  4)  readily 
illustrates  the  scale  of  the  dim  nal  heating  taking  place. 
.A  huge  area  of  SST  warming  ranges  I  rom  0.0  to  0  x  ( 
and  covers  a  section  surrounding  tlie  /one  of  ma.ximum 
warming.  A  light-blue  color  was  purposeU  chosen  to 
offsei  it  from  the  blaek  and  dark -gray  shades  used  for 
cloud-contaminated  pixels.  This  coloration  permits  the 
viewer  easy  visual  access  to  the  image  content.  (Note 
that  there  is  little  fading  from  one  color  to  another  in 
the  color  wedge,  especially  between  1  and  6'C, 
required  to  enhance  each  1  .()°C  day/  night  temperature 
warming  bin). 

Close  examination  is  needed  to  spot  diurnal  heating 
from  1.0°  to  2.0°C,  even  though  it  has  been  colored 
dark-blue  to  offset  it  from  the  previous  light -blue  eolor 
and  the  following  dark-green  (2.0  to  3.0' C)  shade. 
This  1  C  bin  is  scattered  about  the  itnage  but  is  usually 
concentrated  near  the  maximum  heating  sections. 

The  3.0°  to  5.0°C  diurnal  heating  bins  (yellow  and 
orange)  are  concentrated  in  limited  sections  and  nicely 
outline  regions  of  even  higher  heating.  These  values 
typically  fall  between  25  and  35  ’N  and  cover  virtually 
the  whole  longitudinal  range  This  example  of  diurnal 
heating  is  impressive,  even  without  the  larger  values 
found  within  the  yellow-  and  orange-eolored  pixels. 
Note  that  conservative  cloud  tests  were  applied  to 
ensure  that  few  clouds  contaminated  the  day  night 
images.  Thus,  some  of  the  dark  gray  pixels  located  near 
32°N  and  23°W  may  actually  be  cloud-free  oi  partially 
cloudy,  but  were  masked  out  so  that  a  simple  test  could 
be  applied  to  the  entire  image. 

The  DOSL  contours  (valid  for  the  afternoon  of  the 
loth  and  overlain  in  white  in  Fig.  4)  outline  the  general 
diurnal  heating  pattern.  The  very  abriqit  zone  of 
heating  just  off  the  coast  is  well  covered  with  a  DOSl. 
contour  in  excess  of  1.0 'C.  The  model  toreeast 
correctly  predicts  warming  to  occur  around  25' N 
across  a  wide  ocean  expanse.  1  wo  large  pockets  of 
2.5°C  warming  occur  at  26'N,  2*^  W’  and  at  32°N. 
35°W  in  the  DOSL  product,  where  the  wexternmost 
area  agrees  well  with  the  day  'night  M(  SSI  imagery 
The  eastern  maximum  is  slightly  misplaced  to  the 
southeast. 


Figure  3.  MCSST  image  at  IS37Z  16  August  1986  reveals  extremely  hot  SST  values  off  NW  African  coast.  Clouds  have  been 
blacked  out  using  a  multichannel  cloud-screening  technique.  SSTs  have  been  falsely  colored  according  to  temperature  scale 
in  upper-right  corner.  Note  patch  in  excess  of  26°C  near  3l°N,  16°W. 
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Figure  4.  The  16  August  1986  day-night  MCSST  difference  image,  color  coded  according  to  temperature  scale  in  0.5°C  and 
I.O°C  increments.  Light  blue  represents  diurnal  increases  near  1°C;  yellow  is  near  3  °C.  White  contours  are  FNOC  DOSL 
SST  change  forecast  values  at  0.2°C  intervals.  Note  how  well  elongated  region  off  NW  Africa  correlates  with  DOSL  local 
maximum. 
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The  trend  exhibited  in  Figure  4  represents  excellent 
agreement  between  the  DOST  output  and  the  day/night 
MCSST  imagery  when  considering  the  heating  pattern. 
This  agreement  indicates  that  the  wind  stress  and  the 
surface  heat  flu.xes  are  being  correctly  forecast  by 
NOGAPS  in  terms  of  the  overall  meteorological 
conditions  and  are  incorporated  correctly  by  the  DOST 
model’s  physics.  The  few  regions  of  very  intense 
heating,  dramatically  shown  in  imagery,  reveal  that  the 
2.5°C  upper  bound  within  the  DOST  model  can  be  a 
limiting  factor  in  certain  cases.  This  threshold  exists, 
since  the  model  physics  begin  to  break  down  for  values 
above  2.5''C. 

One  should  also  note  the  grid-point  resolution 
involved  in  the  input  fields.  NOGAPS  has  2.5°  spacing 
and  is  extremely  mismatched  when  confronted  with  the 
realities  of  1-km  or  4-km  imagery.  Thus,  the  smoothing 
of  heating  effects  is  to  be  expected  and  will  continue 
to  be  the  norm  until  the  atmospheric  forcing  and  the 
DOST  model  grid  resolution  are  significantly  reduced. 
In  this  light,  the  model  appears  to  do  quite  well  by 
mapping  out  the  prime  diurnal  heating  spots. 

B.  17  August  1986 

The  following  day,  17  August,  also  exhibited  similar 
diurnal  warming  features,  but  not  with  quite  the 
widespread  intensity  evident  on  August  16.  Figure  5 
shows  a  large  east-west  section  of  DOSL-predicted 
afternoon  warming  that  closely  agrees  with  the 
0.0°-0.5°C  (light  blue)  MCSST  day/night  difference. 
The  eastern  segment  is  not  quite  as  sharply  defined  as 
in  the  previous  day,  but  some  cloud-free  pixels  have 
likely  been  masked  out  due  to  high  surface  reflectivity 
values  caused  by  sunglint  near  30°N,  15°W. 

The  DOST  maximums  near  30°N,  18°W  and  33°N, 
32°W  are  somewhat  cloudy;  thus,  good  comparisons 
are  difficult.  However,  the  imagery  clearly  indicates 
several  zones  of  3°  and  4°C  heating  between  the  two 
westernmost  DOST  warming  areas.  Thus,  the  overall 
trend  appears  to  be  correct,  but  individual  locations 
are  hard  to  verify  due  to  spotty  cloud  contamination. 

It  should  be  noted  that  although  the  DOST  model 
has  a  2.5°C  threshold,  the  day/night  MCSST 
maximum  values  of  3°-4°C  compare  well  with  the 
earlier  results  of  Cornillon  and  Stramma  (1985)  in 
the  Sargasso  Sea.  Both  areas  are  affected  by  the  fringe 
of  the  Azores-Bermuda  high,  but  on  opposite  sides. 
The  eastern  North  Atlantic  has  the  advantage  that  no 
major  current  system  with  large  SST  signals  runs 
through  the  area.  The  Gulf  Stream  forces  investigators 
to  move  to  the  south,  away  from  the  disturbing  changes 
that  can  occur  due  to  the  energetic  Gulf  Stream  cold- 
core  rings. 

C.  6  August  1987 

This  example  involves  a  situation  in  the  central  North 
Atlantic  using  NORDA-acquired,  4-km  GAC  data. 


Figure  6  shows  the  data  at  about  4  km  resolution 
(considerably  better  than  previous  cases)  and  permits 
us  to  focus  on  a  zone  within  the  Azores-Bermuda  high. 
A  pronounced  region  of  0.0'’-0.5°C  warming  is  out¬ 
lined  in  light  blue  and  agrees  well  with  DCiSl  cont(,urs 
in  shape  and  extent.  More  imi)oitantly,  the  main  area 
of  inten.se  heating  falls  well  within  the  2.U"'-2.5°C 
DOST  values. 

Diurnal  heating  values  ot  3’  4T.  align  along 
29°-30°N  and  30°-45°W.  These  nurnbeis  exceed  the 
DOST  threshold  (as  in  the  previous  cases),  but  speak 
well  for  the  model,  coiisideiing  the  resolution  for  input 
"ling  functions.  The  fact  that  clouds  lie  to  cither  side 
ihe  maximum  warming  indicates  that  NOGAPS 
appears  to  have  done  well  in  torecasting  the 
meteorological  conditions.  The  extra  region  of  high 
values  in  Figure  6  suggests  that  NOG.APS  is  in  error. 
A  band  of  clouds  oriented  northeast-southwest 
separates  the  two  high-value  sections,  but  DOST  misses 
the  northwestern  segment  by  labeling  it  as  0.2°-0.5°C 
instead  of  the  2°-4°C  heating  present  in  the  imagery. 

The  corresponding  surface  winds  at  OOZ  6  August 
1987  are  outlined  in  the  FNOC  Marine  Wind  Analysis 
(Fig.  7).  This  analysis  is  separate  from  the  NOGAPS 
fields  and  includes  a  combination  of  ship-,  island-, 
and  satellite-derived  wind  retrievals.  Isobars  at 
4-millibar  (mb)  contour  intervals  are  superimposed 
with  high-  (H)  and  low-  (L)  pressure  centers  and  with 
vectors  (half  barb  =  5  kt,  full  barb  =  10  kt,  etc.). 

As  expected,  the  Bermuda  High  dominates  the 
meteorological  pattern  throughout  most  of  the  North 
Atlantic  from  25°N-35°N.  Winds  near  15-20  kt 
surround  the  region,  but  much  calmer  conditions  reside 
in  the  high-pressure  center  as  warm,  dry  air  subsides 
from  above.  This  situation  helps  to  reduce  tropical 
cloud  development  and  promotes  diurnal  heating. 

These  light  winds  are  also  verified  by  the  Special 
Sensor  Microwave/Imager  (SSM/I)  launched  June 
1987  aboard  a  Defense  Meteorological  Satellite 
Program  spacecraft.  This  passive  microwave  sensor  is 
similar  in  frequency  and  purpose  to  the  Nimbus-7 
scanning  multichannel  microwave  radiometer,  but  has 
finer  spatial  resolution  (25  km)  and  better  on-board 
calibration.  It  can  measure  surface  wind  speed  (2  m/s), 
columnar  water  vapor  (0.3  g/cm’),  and  sea-ice 
concentration  (10  percent),  among  other  parameters 
(Hollinger,  1989). 

The  SSM/I  can  measure  surface  winds  in  essentially 
all  weather  conditions.  Heavy  rain  cells  can  cause 
erroneous  retrievals,  but  these  values  can  be  detected 
by  flagging  regions  of  heavy  precipitation,  another 
parameter  sensed  by  the  SSM/ 1.  The  l394-km  SSM/I 
swath  permits  wide  sector  viewing  during  its  14  polar 
orbits.  However,  data  gaps  exist  below  50°,  since 
consecutive  passes  do  not  overlap.  The  gap  increases 
to  1000  km  at  the  equator  and  is  tlicn  filled  in  largely 
by  orbital  data  12  hours  later,  ('omplete  global 
coverage  occurs  every  36  hours. 
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Figure  5.  The  17  A  ugusi  1986  day-night  MCSST difference  image,  color  coded  according  to  temperature  scale.  Genera!  agreement 
in  western  portion,  but  DOSL  values  are  too  high  in  area  of  previous  maximum,  as  seen  in  Fig.  4. 
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Figure  6.  The  6  August  1987  day-night  MCSST  difference  image,  color  coded  according  to  temperature  scale.  This  image  focuses 
on  area  near  center  of  subtropical  high  to  reveal  intense  diurnal  heating  in  excess  of  4°C,  but  is  well  mapped  by  the  DOSL 
contours.  Note  this  figure  uses  higher  spatial  resolution  than  incorporated  in  Figs.  4  and  5. 
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Figure  7.  NO  A  A  National  Hurricane  Center  surface  weather  analysis  for  1200Z  16  August  1986.  Pressure  contours  at  4-rnb 
intervals;  wind  vectors  of  one-half  barb  =  5  kt  2  m/s)  and  one  barb  =  10  kt  (5  m/s).  The  two  sets  of  numbers  next  to 
each  ship  weather  observation  are  air  temperature  (°C)  and  surface  pressure  (multiply  value  by  0.1  and  add  1000  mb  to  get 
correct  value).  Sky  cover  is  marked  by  the  black  portion  of  observation  circle  (i.e.,  open  circle  with  a  line  denotes  no  clouds, 
half  a  black  circle  indicates  presence  of  50Vo  dear  skies).  Note  dear  skies  in  center  of  subtropical  high  near  30°N,  40°  fV. 


LONGITUDE 


Figure  8.  SSM /! surface  winds  on  6  August  1987 for  0700Z  to 
lOOOZ.  The  25-km  wind  retrievals  have  been  mapped  into  a  two- 
dimensional  field  to  examine  the  general  wind  speed  patterns. 
A  large  zone  of  winds  less  than  10  kt  (5  m/s)  parallels  30°N. 


Figure  8  represents  several  passes  covering  the  DOSL 
image  domain  for  5  August  1987  (no  SSM/I  values 
available  on  6  August  due  to  satellite  coverage  gaps). 
The  wind  speed  retrievals  are  contoured  at  5-m/s 
intervals  and  are  shaded  to  help  the  viewer  easily 
recognize  the  zones  of  constant  wind  speed.  Light 
winds  dominate  the  region  and  provide  accurate 
backup  information  to  the  marine  wind  analysis. 

VII.  Future  Considerations 

The  NOGAPS  model  has  undergone  considerable 
revision  since  1987  and  now  contains  an  entirely 
new  treatment  for  the  atmospheric  boundary  layer. 
The  new  approach  has  proven  to  be  far  superior  to 
the  earlier  version,  so  heat  fluxes  for  DOSL  input  have 
been  upgraded  significantly  since  this  study  began. 

The  NOGAPS  and  DOSL  spatial  resolution  of 
2.5°  is  quite  coarse  when  attempting  to  nail  down 
diurnal  heating  features  that  can  be  very  localized. 
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Future  reductions  in  model  grid  spacing  will  help  to 
solve  this  problem. 

Sampling  day  and  night  SSTs  to  validate  DOSL 
results  with  one  polar  orbiter  is  difficult.  Any  follow-on 
studies  should  seriously  consider  the  use  of 
geostationary  data  that  provide  views  every  30  minutes 
and  that  remedy  some  of  the  cloud  contamination 
problems  caused  by  compositing  over  short  intervals 
as  clouds  move.  The  critical  element  is  the  accuracy 
of  SSTs  from  the  Geostationary  Operational 
Environmental  Satellites  (GOES).  The  VISSR 
Atmospheric  Sounder  has  reached  accuracies  of  less 
than  1°C,  but  improvements  are  needed  before  it  is 
equivalent  to  MCSSTs.  This  will  likely  come  about  with 
GOES-NEXT  in  1990-91,  where  multi-infrared 
channels  and  10-bit  digitization  will  open  up  more 
applications  for  these  platforms. 

VIII.  Summary  and  Conclusions 

As  shown  by  the  work  of  Cornillon  and  Stramma 
(1985),  Stramma  et  al.  (1986),  and  many  others,  diurnal 
SST  variability  is  ubiquitous  and  widespread  during 
the  spring  and  summer.  Forced  by  day/night  variations 
in  the  net  surface  heat  flux  and  strongly  modulated 
by  the  magnitude  of  the  surface  wind  stress,  this 
phenomenon  dominates  short-term  variability  in  upper- 
ocean  thermal  stratification  (Price  et  al.,  1986),  and 
thus  critically  impacts  on  surface-duct  acoustics. 

Diurnal  SST  variability  is  also  important  when  using 
and  interpreting  satellite  SST  data.  Satellites  measure 
the  skin  temperature  of  the  ocean,  which  may  or  may 
not  be  characteristic  of  the  average  temperature  above 
the  seasonal  thermocline,  depending  on  the  amplitude 
and  phase  of  the  diurnal  SST  cycle  at  the  time  of 
measurement.  Furthermore,  because  the  diurnal  SST 
response  is  such  a  strongly  nonlinear  function  of 
the  wind  speed,  smoothly  varying  patterns  of  the 
synoptic  wind  field  tend  to  produce  sharp  horizontal 
gradients  in  diurnal  SST  features,  as  demonstrated  by 
the  satellite  imagery  presented  here.  These  sharp 
horizontal  gradients  would  be  easy  to  misinterpret  as 
water  mass  boundaries. 

Price  et  al.  (1986)  developed  scaling  relationships 
relating  the  amplitude  of  the  upper  ocean’s  diurnal 
response  to  wind-stress  and  heat-flux  parameters.  They 
verified  these  relationships,  referred  to  as  the  DOSL 
formulation,  by  using  observed  winds,  heat  fluxes,  and 
SSTs  at  a  moored  buoy. 

Results  from  a  real-time  synoptic  implementation 
of  the  DOSL  model  at  FNOC  were  presented  and  the 
predictions  of  this  model  were  verified  with  satellite 
MCSST  data.  The  case  studies  discussed  here  were 
biased  to  the  extremely  warm  diurnal  heating  events, 
since  the  infrared  imagery  must  obviously  be  cloud- 
free  for  use  and,  thus,  large-magnitude  diurnal  heating 
has  an  increa.sed  chance  to  occur.  The  DOSL  correctly 


defined  the  shape  or  pattern  of  each  heated  region,  but 
the  upper  bound  on  the  warming  is  too  low  by  virtue 
of  the  self-imposed  upper  bound  of  2.5°C.  Day/night 
MCSST  differences  reached  3°-5°C  in  several  cases 
when  light  winds  combined  with  intense  summertime 
insolation  in  the  subtropical  highs. 

IX.  Recommendation 

DOSL  has  proven  to  handle  these  extreme  cases  well 
in  a  qualitative  manner  and  thus  can  be  used  for  several 
immediate  applications.  This  report  recommends  that 
the  DOSL  model  be  declared  operational  and  used  to 
fill  the  gap  existing  in  forecasting  daily  SST  diurnal 
warming  events. 
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